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We develop one-loop effective vertices and propagators in the linear sigma model at finite 
temperature satisfying the chiral Ward identities. We use these in turn to compute the pion 
dispersion relation in a pion medium for small momentum and temperatures on the order of 
the pion mass at next to leading order in the parameter m^./47r 2 /^ and to zeroth order in 
the parameter (m n /m,j) 2 . We show that this expansion reproduces the result obtained from 
chiral perturbation theory at leading order. The main effect is a perturbative, temperature- 
dependent increase of the pion mass. We find no qualitative changes in the pion dispersion 
curve with respect to the leading order behavior in this kinematical regime. 

PACS numbers: ll.10.Wx, 11.30. Rd, 11.55.Fv, 25.75.-q 

The propagation properties of pions in a thermal hadronic environment are a key ingredient for the proper 
understanding of several physical processes taking place in dense and hot plasmas. The scenarios include 
dense astrophysical objects such as neutron stars -where the inclusion of mesonic degrees of freedom is essen- 
tial to determine the equation of state of matter in the inner core- and the evolution of the highly interacting 
region formed in the aftermath of a high-energy heavy-ion collision. In the latter context, knowledge of the 
pion dispersion relation would help to properly account for the component of the dilepton spectrum origi- 
nated in the hadronic phase, below the critical temperatures and densities for the formation of a quark-gluon 
plasma (QGP) Q, and to shed light on possible collective surface phenomena which are speculated to occur 
as a consequence of the mainly attractive interaction between pions and nuclei |^] . 

In order to account for the hadronic degrees of freedom at temperatures and densities below the QCD 
phase transition to a QGP, one needs to resort to effective chiral theories whose basic ingredient is the 
fact that pions are Goldston bosons, originated in the spontaneous breakdown of chiral symmetry. Chiral 
perturbation theory (xPT) is one of such effective theories that has been successfully employed to show the 
well known result that at leading perturbative order and at low momentum, the modification of the pion 
dispersion curve in a pion medium is just a constant, temperature dependent, increase of the pion mass ||. 
Any deviation from this leading order behavior in perturbation theory can only be looked for at higher 
orders. Second order calculations of the pion dispersion curve can be found in Refs. 

The simplest realization of chiral symmetry is nevertheless provided by the much studied linear sigma 
model which possesses the convenient feature of being a renormalizable field theory, both at zero [|| and 
at finite temperature M. This model has lately received a boost of attention in view of recent theoretical 



results H and analyses of data Q that seem to confirm, though not without controversy |f(J , that a broad 
scalar resonance, with a mass in the vicinity of 600 MeV -that can be identified with the cr-meson- indeed 
exists. 

In addition to the above mentioned characteristics, the linear sigma model also reproduces -as we will 
show- the leading order modification to the pion mass in a thermal pion medium (for temperatures on the 
order of the pion mass), when use is made of a systematic expansion in the ratio (m^/iriu) 2 at zeroth order, 
where m CT , m n are the vacuum sigma and pion masses, respectively. This important point, often missed in 
the literature flf|| , provides the connection between the linear sigma model and xPT at finite temperature. 

In this work, we use the linear sigma model to construct effective vertices and propagators to one loop. 
We use these in turn to compute the one loop modification to the pion propagator in a pion medium. We 
require that the effective vertices and propagators thus constructed satisfy the chiral Ward identities [Q. 
The net result is a next to leading order calculation in the parameter m 2 /47r 2 / 2 , where is the pion decay 
constant and to zeroth order in the parameter (m v /m a ) 2 . From the effective pion propagator, we explore 
the behavior of the pion dispersion curve at small momentum and for temperatures on the order of the 
pion mass. Though the momentum dependence of the correction term is non-trivial, we find no qualitative 
difference from the leading order result in this regime. Possible effects due to a high nucleonic density as 
well as general details of the calculation are left out for a future work. 

The Lagrangian for the linear sigma model, including only the mesonic degrees of freedom and after the 
explicit inclusion of the chiral symmetry breaking term, can be expressed as |(| 

C =l [{dv? + {do? - mW rnW] - >?fM° 2 + **) ~ + tt 2 ) 2 , (1) 

where tt and a are the pion and sigma fields, respectively, and the coupling A 2 is given by 

From the Lagrangian in Eq. (Q), one obtains the Green's functions and the Feynman rules to be used in 
perturbative calculations, in the usual manner. Thus, in particular, the bare pion and sigma propagators 
A ff (P), A CT (Q) and the bare one-sigma two-pion and four-pion vertices r^ 2 , 1^ are given by (hereafter, 
capital letters are used to denote four momenta) 

i^(P)S ij = p2 % 2 5 ij 
r — m£ 

iT W = _ 2i x 2 (S»S M + S ik 5 jl + 5 tl S lk ) . (3) 

These Green's functions are sufficient to construct the modification to the pion propagator, both at zero and 
finite temperature, at any given perturbative order. 

An alternative approach consists on exploiting the relations that chiral symmetry imposes among different 
n-point Green's functions. These relations, better known as chiral Ward identities (xWIs), are a direct 
consequence of the fact that the divergence of the axial current may be used as an interpolating field for 
the pion. For example, two of the xWIs satisfied -order by order in perturbation theory- by the functions 
A n (P), A a (Q), Tg and T^f are @ 
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/*r« w (; 0, P 1: P 2 , P 3 ) = (P 1 ;P 2 , P 3 )S^ + T l ^(P 2 ;P 3 , P 1 )S ik + r{ k (P 3 ;P 1 ,P 2 )S il 

UT%(Q;0,P)= [A- 1 (Q)~A- 1 (P)]S t \ (4) 

where momentum conservation at the vertices is implied, that is Pi + P 2 + P 3 = and Q + P — 0. Therefore, 
any perturbative modification of one of these functions introduces modifications in other, when the former 
are related to the latter through xWIs. 

At one loop and after renormalization, the sigma propagator is modified by finite terms. At zero tempera- 
ture this modification is purely imaginary and its physical origin is that a sigma particle, with a mass larger 
than twice the mass of the pion, is unstable and has a (large) non-vanishing width coming from its decay 
channel into two pions. At finite temperature the modification results on real and imaginary parts. The real 
part modifies the sigma dispersion curve whereas the imaginary part represents a temperature dependent 
contribution to the sigma width. 




FIG. 1. Dominant contribution to the sigma self-energy at one loop. 

After renormalization, the leading one-loop contribution to the sigma self-energy, in an expansion in the 
parameter {m^ /m a ) 2 , can be shown Jl5| ] to arise from the diagram depicted in Fig. 1. In the imaginary-time 
formalism of Thermal Field Theory (TFT), the expression for this diagram is given as 6A 4 /^J(w, q), where 
the function I is defined by 

J(<0. <l)= T Ylj ( 2?r )3 Jl n + k 2 + m 2 ( Wb _ w )2 + ( k _ q )2 + m 2 ' ( 5 ) 

Here u> = 2mirT and oj n — 2nirT (m, n integers) are discrete boson frequencies and q = |q|. From Eq. (^|) 
we obtain the time-ordered version P of the function /, after analytical continuation to Minkowski space. 
The imaginary part of P is given by 

Im/*(go, q) = -> Qo + ie, q) - I(icJ q a - ie, q)] 

At 

1 f IT /l _ p-w+(9o.?)/ r \ 1 



16tt \ y ' q \1 - e -u-(io,i)/T 

where Q 2 = q^ — q 2 , e and O are the sign and step functions, respectively, and the functions a and ui± are 
given by 



a(Q 2 ) 



Am 2 



w± («),«) = = (7) 



whereas the real part of P at Q — is given by 
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1 r°° dk 

ReP(0) = --l -[1 + 2/(^)1, (8) 



where Ek = \/k 2 + m 2 and the function / is the Bose-Einstcin distribution 



Therefore, the one-loop effective sigma propagator becomes 

>K(Q) = Q2 _ mi+ \ ximt{Q y do) 

The temperature-independent infinities are absorbed into the redefinition of the physical masses and coupling 
constants by the introduction of suitable countertcrms in the usual manner. 

In order to preserve the xWIs expressed in Eq. (||) , the corresponding one- loop effective one-sigma two-pion 
and four-pion vertices become 

iT* 2 ij (Q;P 1> P 2 ) = -2i\ 2 f^ [1 - 3A 2 /*(Q)] 
iT *W(. Pl P ^ p 3j Pi) = _ 2lX 2 1 [ X _ 2,\ 2 I\P 1 + P 2 )] S^5 kl 

+ [l-3A 2 /*(Pi+P 3 )] S ik S jl 

+ [i - i\ 2 i t {p 1 + p^] s u s ik } . (n) 

It can be shown that the functions in Eq. (|ll| ) arise from considering all of the possible one-loop 
contributions to the one-sigma two-pion and four-pion vertices, when maintaining only the zeroth order 
terms in a systematic expansion in the parameter (m^ / m a ) 2 . 
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(a) (b) (c) 

FIG. 2. Dominant contributions to the pion self-energy at one loop in the effective expansion. The heavy dots 
denote the effective vertices and propagators. 

We now use the above effective vertices and propagator to construct the one-loop modification to the pion 
self-energy. The leading contributions come from the diagrams depicted in Fig. 2. These arc, explicitly 



U a (P) *« = \ 2 6 lj jf * * 4 R2 z m 2 {5 - 9J*(0) - 6X 2 I l (P + K)} 





2xij (Q\ 2 f%\ f d±K I I [l-2,X 2 I\Q)f 



4A 2 /2\ r d 4 K 1 I [1-3X 2 P(P + K)f 



n.(P) S - - -X S' 3 I ) ^ J^J^ - ^ ; - - > ■ (12) 



7T- 
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where the subindex /3 means that the integrals are to be computed at finite temperature. We now expand 
the denominators in the second and third of Eqs. (p"2|), keeping only the leading order contribution when 
considering T and P as small compared to m a . Adding up the above three terms and to zeroth order 
in (m- n jm a )' 1 where 



the pion self-energy can be written as 

,2 



2/r 



d 3 k 1 f fP 2 + K 



(ml\ ^ f d 3 k 1 
( )_ \m) \J (2tt)3 K 2 + m 2 f 

g|) [9/*(0) + 62*(P + JO]} , (14) 

where we carry out the calculation in the imaginary-time formalism of TFT with if = (a> n ,k) and P — (ui,p). 
The pion dispersion relation is thus obtained from the solution to 

P 2 + ml + ReH(P) = , (15) 

after the analytical continuation iu> — ► p + ie. 

As it stands, Eq. (O) contains a temperature-dependent infinity coming from the product of the vacuum 
piece in Re/ (0) and the temperature-dependent piece in the indicated integral, as well as vacuum infinities. 
However, it is well known that finite temperature does not introduce new divergences and that whatever 
regularization and renormalization is needed at zero temperature, it will also be necessary and sufficient at 
finite temperature. Let us recall that at one loop, renormalization is carried out by introducing appropriate 
counterterms into the original Lagrangian and that at this level, no temperature-dependent infinities arise. 
Going to two loops, one will always encounter temperature-dependent infinities in integrals involving only 
the bare terms of the original Lagrangian. The remarkable property of renormalizable theories is that these 
temperature-dependent infinities are exactly canceled by the contribution from the loops computed by using 
the counterterms that were introduced at one loop. The above was explicitly shown for the case of the 
self-energy in the 4> 4 theory by Kislinger and Morley |Q and for the sigma vacuum expectation value in the 
linear sigma model by Mohan JtJ . For this reason, here we omit the renormalization procedure and consider 
only the finite, temperature-dependent terms throughout the rest of the calculation and refer the reader to 
the cited references for details. 

Let us first look at the dispersion relation at leading order. After analytical continuation all the terms are 
real. The integrals involved are 



(2tt) 3 K 2 + ml 2ir 2 J Q E k 
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2tt 2 

d 3 k K 2 (ml 



(2ir) 3 K 2 +ml * \2ir 2 

where g is a dimensionless function of the ratio T/m^ and the arrows indicate only the temperature depen- 
dence of the expressions. Thus, the dispersion relation results from 
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[1 + 2 ^ g(T/m„)] {pi - p 2 ) - [1 + 3 £ g(T/m*)] m 2 = , (17) 

with £ = m 2 /47r 2 / 2 <C 1. For T ~ m w , gr(T/m T ) ~ 1, therefore, at leading order and in the kinematical 
regime that we are considering, Eq. (FTTH can be written as 



pl=p 2 +ml[l+Zg{T/m„)] , (18) 

which coincides with the result obtained from xPT fs|. 

We now look at the next to leading order terms in Eq. fll4| ) . The first of these is purely real and represents 
a constant, second order shift to the pion mass squared 

- • (S§)' r ? / (WiS; ^ 3 e 9(T/ "" ) " (T/ "" ) m - < (19) 

where h is a dimensionless function of the ratio T j defined by 

C°° dk 

h(T/m„) ee J ^rf(E k ) . (20) 



The remaining term in Eq. (14) shows a non-trivial dependence on P. It involves the function S defined by 

The sum is performed by resorting to the spectral representation of /' and (K 2 + to 2 ) -1 . Thus, the real 
part of the retarded version of S, after analytical continuation is 

ReS r (p ,p) = i [S(iu -^p + ie,p) + S{iu -> p - ic,p)] 



2tt e(k' ) 5[k' 2 - (k - p) 2 - to 2 ] 2 Im/*(fc , k) 



(22) 



p -k -k' 

where V represents the principal part of the integral. 

The integration in Eq. ( |22] ) can only be performed numerically. Figure 3 shows a plot of the temperature- 
dependent terms of the function S(po,p) = — (247r 4 )ReS' , '(pojP) f° r T = and po taken as the solution 
of Eq. dig ) up to the P-independent terms of Eq. ( |l4| ) . ^From this figure we see that S is non- vanishing at 
p — and thus it also contributes to the perturbative increase of the pion mass. We also notice that S is 
a monotonically increasing function of p in the kinematical range considered. However, this increase is not 
strong enough to cause a qualitatively significant change on the shape of the dispersion curve, given that 
this function contributes multiplied by the perturbative parameter £ 2 . 
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FIG. 3. S(po,p)/m% as a function of p/m^ for T — m^. The function shows a mild monotonical increase in the 
kinematical region considered. This increase is not enough to change the qualitative behavior of the dispersion curve, 
given that S contributes at second order in the expansion parameter f. 

Including all the terms, the dispersion relation up to next to leading order, for T ~ and in the small 
momentum region is obtained as the solution to 

vl = P 2 + [l + es(TK) + y(?(TK) [9h(T/m w ) ~ 4g{T/m w )]\ml + fSfap) . (23) 





FIG. 4. Pion dispersion relation obtained as the solution to Eq. ( p3| ) for T = (upper curve). Shown is also 
the dispersion relation obtained by ignoring the term S(po,p) (lower curve). 



Figure 4 shows the dispersion relation obtained from Eq. (E3|) for T = m, where we also display the 
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solution without the term £ S(po,p). As mentioned before, inclusion of this last term does not alter the 
shape of the dispersion relation in this kincmatical regime. 

In conclusion, we have shown that in the linear sigma model at finite temperature, the one-loop modifica- 
tion of the sigma propagator induces a modification in the one-sigma two-pion and four-pion vertices in such 
a way as to preserve the %WIs. Strictly speaking, chiral Ward identities fully constrain vertices only in the 
limit where all the external momenta are zero Q. Nevertheless, we have also checked that these functions 
arise from considering all of the possible contributions to the vertices of interest, when maintaining only the 
zeroth order terms in an expansion in the parameter (m ff / 'to ct ) 2 p3[ . We have used these objects to compute 
the next to leading order correction to the pion propagator in a pion medium for small momentum and for 
T ~ m n . We have shown that the linear sigma model yields the same result as xPT at leading order in the 
parameter £ = m 2 /47r 2 / 2 when use is made of a systematic expansion in the parameter (rn n / 'to ct ) 2 at zeroth 
order. This result was to be expected since the kinematical regime we consider is that where the temperature, 
the pion momentum and the pion mass are treated as small quantities, such as in the case of %PT. The main 
modification to the pion dispersion curve in the considered kinematical regime is a perturbative increase of 
the in-medium pion mass. The shape of the curve is not significantly altered. For a sigma meson with a mass 
on the order of 600 MeV and in the mentioned kinematical regime, corrections in powers of the parameter 
(m„/ ma) 2 are small. If a lighter sigma meson exists, corrections of order (m^/rria-) 2 could be important. It 
remains to include possible effects introduced by a high nucleonic density. Finally, it is interesting to note 
that the formalism thus developed could be employed to explore the behavior of the pion dispersion curve 
in the large momentum region where, in principle, the lowest order xPT Lagrangian cannot be used. These 
issues will be treated in a following up work. 

Support for this work has been received in part by CONACyT Mexico under grant numbers I27604-E, 
29273-E and 32279-E. 



[1] C. Gale and J. Kapusta, Phys. Rev. C 35, 2107 (1987). 

[2] E. V. Shuryak, Phys. Rev. D 42, 1764 (1990). 

[3] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987). 

[4] A. Schenk, Phys. Rev. D 47, 5138 (1993); R. D. Pisarski and M. Tytgat, ibid 54, R2989 (1996); D. Toublan ibid 
56, 5629 (1997). 

[5] The two loop computation of the imaginary part of the self energy was carried out by J. L. Goity and H. 
Leutwyler, Phys. Lett. B 228, 517 (1989). 

[6] B. W. Lee, Chiral Dynamics (Gordon and Breach, 1972). 

[7] L. R. R. Mohan, Phys. Rev. D 14, 2670 (1976). 

[8] N. A. Tornqvist and M. Roos, Phys. Rev. Lett. 76, 1575 (1996); ibid 77, 2333 (1996); ibid 78, 1604 (1997). For a 
recent discussion on the status of the sigma meson see for instance What do we know about sigma?, Workshop on 
Hadron Spectroscopy session III, Frascati Italy, March (1999): Frascati Physics Series XV, 75-252 (1999), Edited 



8 



by T. Bressani, A. Feliciello and A. Fillipi. 

[9] M. Svec, Phys. Rev. D 55, 5727 (1997). 

[10] N. Isgur and J. Speth, Phys. Rev. Lett. 77, 2332 (1996); M. Harada, F. Sannino and J. Schechter, ibid 78, 1603 
(1997). 

[11] See for example C. Contreras and M. Loewe, Int. J. Mod. Phys. A 5, 2297 (1990). 

[12] This approach has been previously used to study low energy pion scattering at zero temperature, see J. L. Lucio, 
M. Napsuciale and M. Ruiz-Altaba, The linear sigma model at work: Successful postdictions for pion scattering, 
hep-ph/9903420 . 

[13] Work in progress. 

[14] M. B. Kislinger and P. D. Morley, Phys. Rev. D 13, 2771 (1976). 



9 



